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Abstract 

We compute leading-power fragmentation corrections to J/ij) photoproduction at DESY HERA, 
making use of the nonrelativistic QCD factorization approach. Our calculations include parton 
production cross sections through order a)?, fragmentation functions though order a 2 , and lead¬ 
ing logarithms of the transverse momentum divided by the charm-quark mass to all orders in 
a s . We find that the leading-power fragmentation corrections, beyond those that are included 
through next-to-leading order in a s , are small relative to the fixed-order contributions through 
next-to-leading order in a s . Consequently, an important discrepancy remains between the experi¬ 
mental measurements of the J/U photoproduction cross section and predictions that make use of 
nonrelativistic-QCD long-distance matrix elements that are extracted from the J/fj hadroproduc- 
tion cross-section and polarization data. 

PACS numbers: 12.38.Bx, 12.39.St, 14.40.Pq, 13.60.-r 
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I. INTRODUCTION 


Q, 


According to the nonrelativistic QCD (NRQCD) factorization conjecture [lj, the inclusive 
production cross section to produce a quarkonium H in a collision of particles A and B can 
be written as 

d<JA+B^H+X = daA+B^QQ(n)+X {@ H ( n ) ); ( 1 ) 

n 

where the da A+B ^QQ^ +x are the short-distance coefficients (SDCs), which are essentially 
the production cross section of the heavy quark-antiquark pair QQ{n ) in a specific color and 
angular momentum state n , and the (O h (n)) are the corresponding NRQCD long-distance 
matrix elements (LDMEs), which account for the evolution of the QQ(n) pair into the 
quarkonium. 

The LDMEs have a known scaling with v, the heavy-quark velocity in the quarkonium 
rest frame 1]. Usually, in heavy-quarkonium production phenomenology, v is considered to 
be a small parameter ( v 2 « 0.3 for the J/ip), and the sum over n in Eq. (JT]) is truncated at 
relative order v 4 . For H = J/i^, the sum over n, truncated at order v 4 , involves four QQ 
LDMEs: (O^QS?)), (C^Qsf 1 )), {O j MS l £ ] )), and (0 J ^{ 3 P [ f ] )), where the expres¬ 
sions in parentheses give the color state of the QQ pair (singlet or octet) and the angular- 
momentum state in spectroscopic notation. The color-singlet LDME (C>' 7// ^( 3 S'| 1 ^)) can be 
measured in lattice QCD, can be determined from potential models, or can be computed 
from the J/rp decay rate into lepton pairs. It is not known how to calculate the color-octet 
LDMEs from first principles. As a result, they are usually obtained by comparing Eq. (jT]) 
to measured cross sections. 

SDCs for the J/'if) production cross sections and polarizations have been computed 


through next-to-leading order (N 
tion [2-7] and photoproduction 


, 0 ) 


in the strong-coupling constant a s for hadroproduc- 
10]. The LDMEs that are obtained by fitting the result¬ 


ing cross-section and polarization predictions to the experimental data vary considerably, 


depending on the data that are used in the fits. 


cross-section data 


0 


12] and polarization data 


13 


is possible to fit the hadroproduction 


15] simultaneously [sj], but the resulting 


LDMEs give a prediction for the photoproduction cross section that overshoots the data by 


factors of 4-6 at the highest value o 
cross section has been measured 


the quarkonium transverse momentum px at which the 
. Alternatively, one can extract the LDMEs by compar¬ 


ing the NLO predictions for the hadroproduction and photoproduction cross sections with 
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the experimental data jg], but the resulting LDMEs lead to predictions of large transverse 
polarization in hadroproduction at large pr that are at odds with the experimental data |6]. 

Recently, it was found that fragmentation contributions at leading power (LP) in p t that 
go beyond NLO in a s are important in ,J / iJj hadroproduction [l^. With the inclusion of these 
contributions it was possible, for the first time, to use the LDMEs that have been extracted 
from the hadroproduction cross sections alone to obtain predictions for the J/tp polarization 


at large pr that are near zero and are in agreement with the experimental data 


17]. However, 


these LDMEs, when combined with the NLO SDCs for the photoproduction cross section, 




lead to a prediction that also overshoots the HERA data from the HI Collaboration 
by about a factor of 8 in the highest pt bin for which the cross section has been measured. 

Motivated by this discrepancy between theory and experiment and by the large LP frag¬ 
mentation contributions to J/-0 hadroproduction beyond NLO in a s , we compute in this 
paper LP fragmentation contributions to J/'ijj photoproduction. Our approach is based on 
the method that was described in Ref. jp) ]. 

The remainder of this paper is organized as follows. In Sec. HH we briefly review the 
general method for computing the LP fragmentation contributions that was given in Ref. fr?| . 
Section 1III1 contains a description of the details of the computation of the LP fragmentation 
contribution to photoproduction. In Sec. I1VI we present our numerical results. Finally, in 
Sec.® we summarize and discuss our results. 


II. LP FRAGMENTATION 


At large transverse momentum pT, the contribution at LP in pr to the cross section to 


produce a QQ pair in two-body collisions is given by 


2d, 


21 ] 


^ a A+B^QQ(n)+X ~ d&A+B-ti+X ® A-s.QQ(n) , (2) 

i 

where the daA+B^-i+x are the inclusive parton production cross sections (PPCSs) to pro¬ 
duce a parton i, and the D^qq^ are the fragmentation functions (FFs) for a parton i to 
fragment into a QQ pair in color and angular-momentum state n. At the parton level, be¬ 
fore convolution with parton distribution functions, the LP contribution to dafdp\ depends 
asymptotically on pt as 1 jp\. 

In this paper, we consider gluon fragmentation and light-quark fragmentation. The gluon 
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FFs D g ^QQ^ in Eq. (J2]) are given for the 1 S'o 8 ' channel at order a? (LO) in Refs, 
for the 3 S'f' channel at order a s (LO) in Ref. 0 and at order a 2 S (NLO) in Refs, 
and for the 3 Pf^ channels at order a ^ (LO) in Refs 


23 


22 

25 


24 


, 23], 

261 , 


The light-quark FF D 


q—>QQ(n) 

m the 3 S? channel is given at order al (LO) in Refs. 26-30]. Light-quark fragmentation 
in the other color-octet channels vanishes through order ct^. In the color-singlet channel, 
gluon fragmentation occurs at order a 3 (LO) 31]. Light-quark fragmentation in the color- 
singlet channel vanishes through order a^, but charm-quark fragmentation occurs at order 


a: 


29 


30]. As we will explain in more detail below, an estimate of the size of the color- 


singlet LP fragmentation contribution shows that it is negligible in comparison with the 
overall theoretical uncertainties. Therefore, we will ignore the color-singlet LP fragmentation 
contributions in our numerical analyses. 

The FFs depend on the factorization scale pf. The dependence on pf is gove rned by the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation 
order in a s , the DGLAP equation reads 


32 


35]. At leading 


dtog ^ \D t ) ~ 2n \P„ P ee ) 



( 3 ) 


where D g = D g ^ QQ{ri) , D s = J2f[ D qf ^QQ(n) + D q f ^QQ(n)}, f is the light-quark or light- 
antiquark flavor, the P t j are the splitting functions for the FFs, and rif is the number of 
active light-quark flavors. In order to match what was done in the NLO calculations of 
Refs, [g, 9], we take rif — 3. We resum the leading logarithms of Pt/wIc by solving Eq. (j3J) 
to evolve the FFs from fif — 2 m c to Pf — m T = y/pj> + 4 m%. 

Following Ref. {l^], we combine the LO-plus-NLO SDCs with the LP fragmentation 
contributions according to the formula 


da LP+NLO da 


LP 


dp T 


dpx dpx 


^^NLO _|_ ^^NLO 


dpr 


( 4 ) 


Here, da hP /dpT is the DGLAP-evolved LP fragmentation contribution, da^Lo/dpr is the 
contribution that arises from the LO-plus-NLO SDCs, and da^[ Q /dpT is the contribution 
that is contained in both da LF /dpx and da^Lo/dpr- 
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III. COMPUTATION OF LP SDCS 


In photoproduction at HERA, the incoming electron or positron emits a virtual photon 
that is nearly on mass shell and that subsequently interacts with the incoming proton. There 
are two types of photon-induced processes that contribute to photoproduction cross sections. 
The first is the direct process, in which the virtual photon interacts with a parton in the 
proton electromagnetically. The second is the resolved process, in which the virtual photon 
emits a parton, which then interacts strongly with a parton in the proton. The probability 
for the photon to emit a parton is given by the parton distribution function (PDF) of the 
photon. 

We compute the PPCSs da e+p ^ i+ x for direct and resolved photoproduction to NLO 


accuracy in a s by making use of the EPHOX Fortran code 
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We use the AFG04_BF 


photon PDFs and the CTEQ6M proton PDFs at scale mx- We carry out the computation 
using the same kinematics and cuts as were used by the HI Collaboration in their most 
recent cross-section measurements 19]. The center-of-momentum energy of the ep system 


is i/s = 319 GeV. The cuts on the yp-invariant mass W = \f (p 7 + p p ) 2 and elasticity 
z = Pj/^ ■ p p /p 7 • Pp are given by 60 GeV < W < 240 GeV and 0.3 < z < 0.9. Here, p 7 , 
Pp, and pj are the momenta of the photon, proton, and J/ip, respectively. There is also a 
cut on the invariant mass of the virtual photon Q 2 , which is Q 2 < QL* = 2.5 GeV 2 . The 
photon flux is calculated in EPHOX by making use of the Weizsacker-Williams formula for 
quasireal photons: 


a 


Mx) = ^ 


1 + (1 - x f log QL*(! - x ) 2 (! - x ) 


x ~ mix 2 x 

Here, x = E^/E e , where E, f and E e are the energy of the photon and the electron, respec¬ 
tively; a is the quantum electrodynamics coupling constant; and m e is the electron mass. 


( 5 ) 


IV. NUMERICAL RESULTS 


In Fig. |H we compare our results for da^ 0 /dpx for the direct process with the NLO 
calculation for the direct process. In Fig. [21 we make a similar comparison for the resolved 
process. Here, and throughout the remainder of this paper, we make use of the SDCs through 


NLO that were computed in Refs. [8, 9, 


160 
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FIG. 1: The ratio {dcr^[ 0 / dpp) / {da^LO / dpr) f° r the 1 S'o^, 3 Pf\ and 3 S'[ 8 ^ channels in the direct 
process 7 p —>• J/•ip + X. 
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FIG. 2: The ratio (do^LO / dpr)/ (doNLO/ dpr) for the 1 5q 8 ^, 3 Pf\ and 3 s[ 8 ^ channels in the resolved 
process -yp —> J/ip + X. 
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For the direct process, LP fragmentation has a sizable contribution only for the 3 Sf^ 
channel. In the 1 S'q 8 ' channel, da^ p Q /dpx is very small in comparison to da^Lo/dpr because 
the PPCSs involve only light-quark-initiated processes, which have much smaller partonic 
fluxes than do gluon-initiated processes, and because the FFs contain no distributions that 
emphasize the region near z — 1. In the 3 Pj^ channels, da^ Q /dpx is also small in comparison 
to da^Lo/dpr (less than 1% for p T < 10 GeV) because the PPCSs involve only light-quark- 
initiated processes. In the 3 Pf^ channels, dcr^ p 0 /dpx is opposite in sign to the non-LP part of 
da^Lo/dpr and, as expected, grows in magnitude relative to the non-LP part of da^o/dpx 
as pt increases. However, even at pt = 10 GeV, da^[ 0 /dpx is much smaller than the non-LP 
part of da-^ho/dpT■ Consequently, the ratio (da^ p Q /dpx)/{da^i J o/dpx) becomes increasingly 
negative as pt increases. At very large values of px, at which the magnitude of d.a\^ 0 /d:pr 
becomes comparable to or larger than the magnitude of the non-LP part of cOnlo/ dpr, we 
would expect this ratio to change sign discontinuously and to approach unity. We do not 
show results for the 3 S'{ 1 ^ channel. We defer the discussion of that channel until we discuss 
the sum of the direct and resolved contributions to the cross section. 

For the resolved process, LP fragmentation in the 3 sf^ and 3 pf^ channels approaches the 
NLO calculation as px rises. In the 1 S'q 8 ^ channel, the LP fragmentation contribution is small 
because, unlike the FFs for the 3 S j 8 ^ and 3 Pj 8 ^ channels, the FF for the 1 S'q 8 ^ channel does not 
involve distributions that emphasize the region near z = 1. Again, we defer the discussion 
of the 3 S'{ 1 ^ channel until we discuss the sum of the direct and resolved contributions to the 
cross section. 

In Figs. [3] and U we compare chx LP+NL0 /dpx with da^xo/dpx in each channel for the 
direct and the resolved process, respectively. For the direct process, da LP+NLO /dpx is larger 
than d<j nlo / dpx in the 3 Sf^ channel, while in the 3 pj 8 ^ and 1 S'o^ channels, the differences 
between da LP+NLO /dp T and dcr N Lo/dpx are less than 1% and less than 0.02%, respectively. 
For the resolved process, the difference between da LP+NLO /dpx and do^Lo/dpx is substantial 
in the 3 Sf^ and 3 Pf^ channels. In the 1 S'q 8 ^ channel, dcr LP+NLO /dpx is larger than da^ho/dpx 
by only 7% at px = 7 GeV and by only 10% at px = 10 GeV. 

In Fig. El we compare, for each channel, the sum of the direct and resolved contributions 
to da LP+NLO /dpx with the sum of the direct and resolved contributions to da^ho/dpx- In 
the 3 Sf] channel, dcr LP+NLO /dp T is smaller than da NLO /dp T by 40% for p T < 10 GeV. On 
the other hand, in the 1 and 3 Pj 8 ^ channels, the differences between dcr LP+NLO / dpx and 
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FIG. 3: The ratio (dcr LP+NLO / dpp) / (da^LO / dpr) for the 1 S’o^, 3 Pf\ and 3 sf^ channels in the 
direct process "yp —> J/ip + X. 
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FIG. 4: The ratio (dcr LP+NLO /dpp) / (da^LO /dpr) for the 1 5’q 8 ^, 3 Pf\ and 3 s| 8 ^ channels in the 
resolved process 'yp —>• J / ip + X. 
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FIG. 5: The ratio (dcr LP+NLO / dpx) /(dcNLO/ dpr) for the 1 5g 8 ^, 3 Pf\ and 3 sf^ channels in the sum 
of the direct and resolved processes yp —>• J/i^ + X. 


d&Nho/dpT are negligible. The reason for this is that, owing to the experimental kinematic 
constraints on W and z, the direct process dominates over the resolved process in these 
channels, and the differences between dcr LP+NLO / dp T and da^ho/dpr are small in direct 
production in these channels [ff|. We do not show the contribution of the 3 S'[ 1 ^ channel. We 
have made a rough estimate of the size of the LP fragmentation correction in that channel 
by making use of the gluon and charm-quark FFs at LO in a s . We estimate that, in the 
3 S[^ channel, the LP fragmentation contributions for 9 GeV< pr < 10 GeV are less than 
5% of the sum of the color-singlet direct and resolved contributions. 

Finally, we can examine the effect of the additional LP fragmentation contributions on 
the complete cross section. As we have seen, the additional LP fragmentation contribution 
is a sizable fraction of the rate only for the 3 Sf^ channel. However, the SDC through NLO 


for the 3 S'j 8 ' channel is much smaller than the SDCs through NLO for the 1 S'o° l and 3 Pj 
channels. Hence, the correction from the additional LP fragmentation contributions to the 
J/ijj photoproduction cross section is very small. For example, if we use the color-octet 


1 o[8] 


,[ 8 ] 


LDM 

data 


is that were obtained in Ref. 


17) by fitting the LP+NLO SDCs to CMS [12] and CDF 




11] and if we use the color-singlet LDME that was obtained in Ref. 40| by making use 
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;0 predictions for the «//?/> differential cross section at HERA compared with the 


of the J/il) leptonic decay-rate data, then we find that the difference between da LP+NLO /dpx 
and da-Nho/dpT for J/^ photoproduction is less than 1% for px between 7 GeV and 10 GeV. If 


we use the LDMEs from Ref. 
is less than 3% for the same px range 


16], then the difference between da LP+NLO /dpx and da^Lo/dpr 


In Fig. IH1 we compare the HI data 18, [19] for the J/ijj photoproduction cross section with 


171 ] and 


the LP+NLO prediction that results from the use of the color-octet LDMEs of Ref. 
the color-singlet LDME of Ref. ^0| . The uncertainty band in the LP+NLO prediction comes 
from the uncertainties in the LDMEs, combined in quadrature with the uncertainties in the 

n 

SDCs that were obtained in Ref. [16|] by varying the factorization scale p,f between 2mx and 
ttit/2. As we have already mentioned, the additional LP fragmentation contributions do 
not significantly change the prediction from that of the NLO calculation. Consequently, the 
LP+NLO prediction overshoots the data by about a factor of 8 at the highest value of px 
at which the photoproduction cross section has been measured. 

In order to suppress next-to-leading-power (NLP) contributions, we have compared the 
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theoretical predictions with data only for values of px that are considerably larger than mj 


In Ref. 


4l|, it was suggested that a criterion px > be used in comparing data with 


theory. In Ref. 17], a criterion px > 10 GeV was used in fitting the predictions to the data. 
The highest value of px at which the photoproduction cross section has been measured falls 
slightly short of these criteria. It is possible that NLP corrections and/or power-suppressed 
violations of NRQCD factorization could account for some of the differences between the 
LP+NLO prediction and the experimental data. However, the data and the prediction do 
not seem to be trending toward each other as px increases. 


The color-octet LDMEs of Ref. 


17] were extracted from the prompt J/if) production data, 


which include feeddown from the XcJ an d ip (2S) states. We would expect the corrections to 
the LDMEs from the removal of the feeddown contributions to prompt hadroproduction to 
be of order —30% 


42], and we would expect the corrections from t 


re inclusion of feeddown 


contributions to prompt photoproduction to be of order +15% 


19]. This still leaves a 


substantial discrepancy between the data and theoretical predictions. 


V. SUMMARY AND DISCUSSION 

In this paper, we have computed additional leading-power (LP) fragmentation contribu¬ 
tions to J/ip photoproduction at HERA that go beyond the existing fixed-order calculations 
through NLO in a s . Our computation has made use of parton production cross sections 
through NLO in a s , as implemented in EPHOX, for both direct and resolved photoproduc¬ 
tion. We have included gluon and light-quark fragmentation processes and have used the 
current state-of-the-art fragmentation functions for the various color-octet QQ channels. For 
gluon fragmentation, the fragmentation function for the 3 S'j 8 ^ channel is available through 
NLO in a s , and the fragmentation functions for the 1 S'q 8 ^ and 3 Pp channels are available at 
LO in ot s . For light-quark fragmentation, only the LO fragmentation function for the 3 Sf^ 
channel is available. In addition to making use of the available fixed-order fragmentation 
functions, we have resummed leading logarithms of p^/m 2 to all orders in a s for the channels 
in which fragmentation functions are available. We have estimated the LP fragmentation 
corrections in the 3 S'P channel by making use of the gluon and charm-quark fragmentation 
functions at LO in a s . 

We find that the additional LP fragmentation contributions are important, relative to the 
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fixed-order contributions through NLO in a s , only in the 3 S'® channel in direct production 
and only in the 3 Pf^ channel in resolved production. However, the fixed-order contributions 
from direct production in the 3 S® channel are themselves small, and the contributions from 
resolved production are small in the 3 Pf^ channel (and the 3 S® channel) for the kinemat¬ 
ics and cuts of the HI cross-section measurement. Furthermore, all of the additional LP 
fragmentation contributions in the 3 S® channel are negligible compared to the theoretical 
uncertainty in the sum of the contributions of all of the channels. Consequently, the ad¬ 
ditional LP fragmentation contributions that we have computed have little effect on the 
cross-section prediction. 

If one predicts the photoproduction cross section by using NRQCD long-distance matrix 
elements that are consistent with both the hadroproduction cross-section and polarization 
data, then there is a sizable discrepancy between theory and experiment at the highest 
values of p t at which the photoproduction cross section has been measured. If NRQCD fac¬ 
torization is correct, then one would expect it to hold at values of p T that are considerably 
larger than rrij/f. As we have mentioned, the highest value of Pt at which the photoproduc¬ 
tion cross section has been measured is not very large in comparison with mj/^. Hence, it 
is possible that NLP corrections and/or power-suppressed violations of NRQCD factoriza¬ 
tion could account for some of the discrepancy between theory and the experimental data. 
However, the shapes of the data and the LP+NLO prediction versus pt do not suggest a 
resolution of the discrepancy at larger values of px- Hence, the discrepancy between theory 
and experiment in photoproduction of the Jj\b seems to challenge the validity of NRQCD 
factorization. 

The calculations in this paper include some, but not all, of the nonlogarithmic LP frag¬ 
mentation contributions at next-to-next-to-leading order (NNLO) in a s (order a 3 ). A com¬ 
plete calculation of the nonlogarithmic LP fragmentation contributions through NNLO in 
a s would require the calculation of additional QCD corrections to fragmentation functions 
and parton production cross sections. In the case of gluon fragmentation, which dominates 
in hadroproduction, a complete NNLO calculation in the 3 S® and 3 P® channels would 
require the use of fragmentation functions through NLO and parton production cross sec¬ 
tions through NLO, and a complete calculation in the 3 S® channel would require the use 
of the fragmentation function through NNLO and parton production cross sections through 
NNLO. In the case of light-quark fragmentation, a complete NNLO calculation in the 3 S® 
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and 3 Pf^ channels would require the use of fragmentation functions through NLO and the 
LO parton production cross sections, and a complete calculation in the 3 S'® channel would 
require the use of the fragmentation function through NLO and parton production cross 
sections through NLO. As we have mentioned, the fragmentation functions are all known 
only at LO, except in the case of gluon fragmentation in the 3 S® channel, in which case 
the fragmentation function is known through NLO. Parton production cross sections are 
publicly available through NLO. 

A complete NNLO calculation of the LP fragmentation contributions might be important 
for hadroproduction cross sections and polarizations, and, hence, could affect the extractions 
of the long-distance matrix elements from the hadroproduction data. However, given the 
small sizes of the additional LP fragmentation contributions that we have found in this paper, 
it seems unlikely that these further LP fragmentation contributions would be important for 
photoproduction. 
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